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Abstract: The impedance matching in metamaterial perfect absorbers 
has been believed to involve and rely on magnetic resonant response, with 
a direct evidence from the anti-parallel directions of surface currents in 
the metal structures. Here we present a different theoretical interpretation 
based on interferences, which shows that the two layers of metal structure 
in metamaterial absorbers are linked only by multiple reflections with 
negligible near-field interactions or magnetic resonances. This is further 
supported by the out-of-phase surface currents derived at the interfaces 
of resonator array and ground plane through multiple reflections and 
superpositions. The theory developed here explains all features observed 
in narrowband metamaterial absorbers and therefore provides a profound 
understanding of the underlying physics. 
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1. Introduction 

The demonstration of metamaterial perfect absorbers [1,2] represents one of the most impor- 
tant applications employing the astonishing properties found in metamaterials [3-5 1. The orig- 
inal idea J2) is that, in metamaterials with simultaneous electrical and magnetic resonances, 
both of the effective permittivity e(co) and permeability fi(co) are highly dispersive and can 
be tailored independently. At certain frequencies the effective impedance, which is defined as 
Z((o) = y/fi(co)/ e(co), matches to the free space impedance Zq, and therefore the reflection 
is minimized. If at the same time the metamaterial is also of high loss which causes very low 
transmission, then near-unity absorption can occur within an ultra thin layer of material. Un- 
der such considerations, in a typical metamaterial absorber, the magnetic resonance from the 
bi-layered metal structure J6l|7) is essential, and each layer of metal structure also provides the 
electrical response. As a direct evidence of the magnetic resonance in metamaterial absorbers, 
the surface currents excited in the two metal layers are found to be anti -parallel QQ. 

While all of the above arguments seem to be plausible and accepted by a majority of re- 
searchers I9 WTT1 . fundamental questions remain. In fact, such a metamaterial absorber can 
be equivalent to a single layer of atoms (or molecules), and it is strongly inhomogeneous 
in the wave propagating direction |12|. So it can hardly be considered as an effective bulk 
medium where the constitutive parameters e(ffl) and ji(co) apply. Additionally, Fabry-Perot 
resonance lfT2Hl4l . conventional transmission line model ifTBI . and cavity resonance lfl6l have 
been proposed to explain metamaterial absorbers or antireflection coatings, which also raise 
questions regarding the existence or involvement of the magnetic resonance in metamaterial 
absorbers. 

In this paper, we start with the multiple reflections interference model [13], and show in a 
typical metamaterial absorber there is negligible near-field interaction or magnetic response 
between the neighboring metal structures. The two layers of metal structures can be decoupled 
and the only link is multiple reflections between them. Through combination of numerical 
simulations and analytical calculations, we further derive the surface currents and reveal that 
the anti-parallel directions are a result of interference and superposition, rather than excited by 
the magnetic component of the incident electromagnetic fields. 



2. Models in numerical simulations 



2.1. Coupled system 

Without losing generality, we use a simple and typical metamaterial absorber structure oper- 
ating at terahertz frequencies as an example. It consists of a cross-resonator array and a ground 
plane separated by d = 10 jJ.m thick polyimide dielectric spacer lfl7HT9l . The unit cell is shown 
in Fig. []I a X where the periodicity is p = 100 jim, the length and width of cross wires are 
/ = 90 jJ.m and w = 10 |im, respectively, and the thickness of metal is t = 0.2 jim. Due to the 
presence of the ground plane, the substrate provides the mechanical support only and the choice 
of material does not affect the performance of the metamaterial absorber. We carry out numer- 




Fig. 1. (a) Schematic structure of the metamaterial absorber unit cell, which repeats in x 
and y directions forming a square array with periodicity of p. (b) Unit cell used to obtain 
the reflection and transmission coefficients at air-spacer interface with the cross-resonator 
array, (c) The solid red curve is the reflection \S\\ \ from the metamaterial absorber shown in 
unit cell (a), and the dashed blue curve is the resonant transmission |52i | through the cross- 
resonator array shown in unit cell (b). (d) Surface currents of the metamaterial absorber at 
the peak absorption frequency. The current at the cross flows to left direction, while it flows 
to right direction at the ground plane. 



ical simulations using CST Microwave Studio 2009 H20I . Perfect electric conductor (PEC) is 
used to simulate the metal, while the ohmic loss is effectively put in the polyimide spacer with 
dielectric constant £ sp acei- = 3.1 x (1 +0.07;). Periodic boundary conditions are applied to the 
side walls parallel to z direction, and the ports are at the front and back surfaces of the unit 
cell. The the simulated electric field reflection |Sn| is shown by the solid red curve in Fig.QJc), 
which reveals a nearly zero reflection at a frequency slightly below 1 THz. This also suggests 
that the incident light is almost completely absorbed by the metamaterial in this narrow fre- 
quency band, according to A(co) = 1 — R(co) — T(co), where R(co) = \Su\ 2 is the reflectance, 
and T(co) = \S2i | 2 = is the zero transmittance due the presence of the ground plane. The sur- 
face currents in the cross-resonator array and the ground plane exhibit anti-parallel directions, 
as shown in Fig.Qld). 

2.2. Decoupled system 

The above model treats the metamaterial absorber as a coupled system, i.e., the possible near- 
field interactions and magnetic resonance have been taken into account between the cross- 
resonator array and the ground plane. However, we may also consider the presence of the metal 
cross-resonator array resulting in an impedance-tuned air-spacer interface ||2TI with dramati- 
cally modified complex reflection and transmission coefficients shown in Fig. [2] The ground 
plane, on the other hand, functions as a perfect reflector with reflection coefficient — 1 . In the 
interference model we decouple the metamaterial absorber into two tuned interfaces, respec- 
tively with cross-resonator array and ground plane located at the two sides of the spacer. They 
are only linked by multiple reflections, as shown in the inset to Fig. [2ja), while any near-field 
interaction or magnetic resonance has been neglected. In this model, we need to know the re- 
flection/transmission coefficients at the air-spacer interface with cross-resonator array, which 
are simulated using the unit cell shown in Fig.QJb), where the ground plane and the substrate 
have been removed from the unit cell shown in Fig. [TJ a )- The cross-resonator array reveals 
a dipole resonance as indicated by the deep transmission dip at about 1.1 THz, shown in the 
dashed blue curve in Fig. [TJ C )- Comparing the two curves, the feature of particular interest 
and worthy of mentioning is that near-zero reflection (or unity absorption) of the metamaterial 
absorber occurs at a frequency different from and lower than the resonance frequency of the 
cross-resonator array. 

3. Interference theory 

As shown in the inset to Fig. [2ja), at the air-spacer interface with cross-resonator array, the 
incident light is partially reflected back to air with a reflection coefficient f\2 = r^e"^ 12 and 
partially transmitted into the spacer with a transmission coefficient t\2 = fi2<?' e ' 2 . The latter 
continues to propagate until it reaches the ground plane, with a complex propagation phase 
p = J3r + ifti = yj £ S pacevk()d , where ko is the f ree space wavenumber, j3, is the propagation 
phase, and /3, represents the absorption in the spacer. After the reflection at the ground plane 
and addition of another propagation phase /3, partial reflection and transmission occur again at 
the air-spacer interface with cross-resonator, with coefficients of ?2i = rzie'™ and t2\ = t2\e' e2i , 
respectively. Similar to the light propagation in a stratified media 1221 . the overall reflection is 
then the superposition of the multiple reflections: 

r=rn-~ r-r, (1) 

1 +r 2 ie ,2 P 



where the first term is the reflection directly from the cross-resonator array, and the second term, 
including the "— " sign, is the reflection resulting from superposition of the multiple reflections 



1.0 
0.8 



Q. 

J 0.4 
< 



0.2 - 



0.0 
270 



_ 180 |- 

ui 
a> 

I 90 

CD 



- 



-90 




_L 



0.6 0.8 1.0 1.2 

Frequency (THz) 



1.4 



Fig. 2. (a) Amplitude and (b) phase of the reflection and transmission coefficients at the 
air-spacer interface with cross-resonator array, obtained by numerical simulations using 
the unit cell shown in Fig. |TJb). Inset: Multiple reflections and interference model of the 
metamaterial absorber, where the cross-resonator array is indicated by the dashed line at 
the air-spacer interface. 



between the cross-resonator array and ground plane. The absorptance is then retrieved through 
A (a) = 1 — |?(fi))| 2 since the transmission is zero. 

3.1. Spacer thickness dependent absorption 

Using the reflection and transmission coefficients shown in Fig. [2 the absorptance is calculated 
in the decoupled metamaterial absorber system and the results are shown in Fig. [3] where we 
have purposely varying the spacer thickness d from 4 jim to 16 [Xm. The spacer thickness 
dependent absorptance in the coupled metamaterial absorber is also simulated directly using 
the unit cell shown in Fig. [TJa) and the results are shown in the inset to Fig. [3] which is in 
excellent agreement with the results calculated in the decoupled metamaterial absorber using 
the interference model. With all of these spacer thicknesses, we observe an absorption peak, 
which increases in the beginning and its frequency red-shifts as the spacer thickness increases. 
It reaches near-unity absorption at 1 THz when the spacer thickness is about 10 jim, and then 
the peak value decreases and its frequency continuously red-shifts when the spacer thickness 
further increases. Further numerical simulations and calculations reveal that, with any chosen 
spacer dielectric material, there is always an optimized spacer thickness where the reflection 
can reach zero, i.e. unity-absorption. 

The excellent agreement shown in Fig. [5] validates the interference model we have pre- 
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Fig. 3. Calculated absorptance in the decoupled metamaterial absorber using the interfer- 
ence model for various spacer thicknesses. Insets: The simulated absorptance when treating 
the whole metamaterial absorber as a coupled system. 

sented above. The physical explanation of a metamaterial absorber is then as follows. The 
multiple reflections in the metamaterial absorber, i.e. the second term in Eq. (HJ, construc- 
tively interfere as evidenced by the fact that near 1 THz the phase change of a round trip is 
2j3r + 021 + 180° « 360°. The superposition of the multiple reflections then destructively inter- 
feres with the direct reflection from the air-spacer interface with cross-resonators, i.e. the first 
term in Eq. ([TJ. With the optimized spacer thickness, these two terms cancel each other out 
resulting in zero reflection; while with other spacer thicknesses the amplitude and phase do not 
match and these two terms only partially cancel each other out, resulting in a reduced absorp- 
tion peak and corresponding frequency shift. Obviously, this explanation does not involve any 
near-field interaction or magnetic resonance between the two metal layers in the metamaterial 
absorber. 

3.2. Anti-parallel surface currents 

The remaining question is then why we observe the anti-parallel directions of the excited 
surface currents as shown in Fig. Hid), which has been used as the evidence of a mag- 
netic resonant response J8) similar to fish-net metamaterials flUE). In order to elucidate this 
mystic observation, by considering the boundary conditions at the interfaces we derive the 
excited surface currents based on interference and superposition. When the incident light 
with electric field ££x is first reflected and transmitted at the air-spacer interface with cross- 
resonator array, according to the boundary condition requirements 11231 . it excites a surface 
current J c = [(1 — ?i2)/Zq — fi2/Z s )]£'£x, where Zq and Z s are impedances of the free space 
and spacer, respectively. The transmitted light continues to propagate and reflected by the 
ground plane, and the latter excites a surface current at the ground plane J g = Ifne'P E\/Z s x. 
The reflection from the ground plane propagates back to the air-spacer interface with cross- 
resonator array, where partial reflection and transmission occur which excites a surface current 
J^. = [(1 — r2i)/Z s — t2\ / Zo]t\2e' 2 P (—\)E l c x. This procedure continues and the overall surface 
currents at the cross-resonator array and ground plane are then superpositions of multiple exci- 



tations: 



r total 



l+r 2 ie a P Z s 



£'x. (3) 



With a spacer thickness of = 10 |im, the phase and magnitude of surface currents in Eqs. 
and (O are shown in Fig.HJa) and its inset, respectively, which reveal a phase difference of 
Ay/ = y/ c — y/g = 192° and comparable surface current magnitude at 1 THz. That is, they are 
almost out-of -phase, directly resulted from the interference rather than the magnetic resonance. 
We also notice that, in the calculated results shown in Fig. SJa), at frequencies below 1 THz 
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Fig. 4. (a) Phase spectra of the excited surface currents at the air-spacer interface with 
cross-resonator array (i// c ) and at the ground plane (i^g), as well as their difference (\)f c — 
i//g), calculated through multiple reflections and superposition based on the interference 
model. The corresponding magnitude of the surface currents are shown in the inset, (b) 
Backward radiation field (reflection) calculated through E r °c J| 
spacer thicknesses. 
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the phase difference is ^180°, while at frequencies above the resonance of the cross-resonator 
array (about 1.1 THz) the phase difference is nearly 360°. This agrees well with the results 
in numerical simulations when treating the whole metamaterial absorber as a coupled system. 
Finally, we use the obtained surface currents to calculate the overall backward radiation field 
(i.e. reflection) by E r Jj. 0tal + Jg° tal e ! ^° for different spacer thicknesses. The results shown in 
Fig. 2|b) reveal a reflection dip, in which the spacer thickness dependence of the reflection 
minimum and its frequency is again consistent with the absorption results shown in Fig. [2] 

4. Conclusion 

We have shown that it is not required in a metamaterial absorber to have simultaneous elec- 
tric and magnetic responses, which have been considered as the foundation of metamaterial 
absorbers. We explicitly demonstrate that the assumed magnetic resonance plays a negligible 
role in the impedance matching of metamaterial absorbers to free space. In contract, it is the 
destructive interference between the direct reflection and the following multiple reflections that 
effectively traps light in the metamaterial absorbers and eventually causes the high absorption. 
Based on such an interference model, we derive the surface currents through multiple excita- 
tions and superposition, which clearly show the anti-parallel directions of surface currents in 
the resonator array and ground plane. This profound understanding of the underlying physics 
will undoubtedly provide valuable guidance in the future developments of more advanced meta- 
material absorber based devices for microwave and photonic applications. 
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